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ABSTRACT 

Recent advances in spectropolarimetry have allowed the reconstruction of stellar coro- 
nal magnetic fields. This uses Zeeman-Doppler magnetograms of the surface magnetic 
field as a lower boundary condition. The ZDI maps, however, suffer from the absence 
of information about the magnetic field over regions of the surface due to the presence 
of dark starspots and portions of the surface out of view due to a tilt in the rotation 
axis. They also suffer from finite resolution which leads to small scale field structures 
being neglected. 

This paper explores the effects of this loss of information on the extrapolated 
coronal fields. For this we use simulated stellar surface magnetic maps for two hypo- 
thetical stars. Using the potential field approximation, the coronal fields and emission 
measures are calculated. This is repeated for the cases of missing information due to, 
(i) starspots, (ii) a large area of the stellar surface out of view, (iii) a finite resolution. 

The largest effect on the magnetic field structure arises when a significant por- 
tion of the stellar surface remains out of view. This changes the nature of the field 
lines that connect to this obscured hemisphere. Nonetheless, the field structure in the 
visible hemisphere is reliably reproduced. Thus the calculation of the locations and sur- 
face filling factors of accretion funnels is reasonably well reproduced for the observed 
hemisphere. The decrease with height of the magnetic pressure, which is important in 
calculating disc truncation radii for accreting stars, is also largely unaffected in the 
equatorial plane. The fraction of surface flux that is open and therefore able to supply 
angular momentum loss in a wind, however, is often overestimated in the presence of 
missing flux. 

The magnitude and rotational modulation of the calculated emission measures is 
consistently decreased by the loss of magnetic flux in dark starspots. For very inactive 
stars, this may make it impossible to recover a magnetic cycle in the coronal emission. 
Finite resolution has little effect on those quantities, such as the emission measure and 
the average coronal electron density, that can currently be observed. 



1 INTRODUCTION 



Since magnetic activity was first detected on stars other than 
the Sun ( Babcock|1947 ), a great deal has bean learned about 
the structure of stellar magnetic fields and their effects on 
the evolution of stars. Using the Zeeman effect on magnet- 
ically sensitive spectral lines, it is possible to measure the 
mean surface fields (using Zeeman splitting/broadening) and 
the mean line-of-sight fields (using circular polarization) for 
unresolved stars. 

Significant progress has been made since it has become 
possible to identify individual structures on the surfaces of 
unresolved stars as they rotate. This is done using a tech- 
nique called Doppler-imaging which allows surface tempera- 
ture and chemical composition maps of rapidly rotating stars 
to be reconstructed from rotationally broadened spectral 
lin es dVogt Penrod|1983| . Dopple r-imaging was extended 
by Semel ( 1989 ) to a technique called Zeeman-Doppler imag- 



ing (ZDI) which uses the Stokes I and circularly polarized 
Stokes V parameters to give reliable surface magnetic maps 
(e.g. Donati & Collier Cameron 1997| 
Morin et al.||2008D . 



Donati et al. 2006 



This technique has been used with great success to 
study the kinds of surface magnetic fields that many stars 
possess. For example, [Donati et al. ( ,2006) produced ZDI 
maps for the massive, slowly rotating star r Sco. It was found 
that T Sco possesses a strong, complex, mostly poloidal field 
that is distributed over all latitudes and constant over the 
period of 1.5 years that the star was observed. These ob- 
servations are inconsistent with current dynamo models for 
massive stars. They concluded that r Sco probably possesses 
a fossil- field. 

For some studies, such as those involving the magnetic 
interaction between a young star and its disc, knowledge of 
the 3-dimensional structure of a star's magnetic field is also 
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required. Making a few simplifying assumptions, and using 
ZDI maps as a lower boundary condition, the structure of 
the coronal fields can be extrapolated. The most common 
technique that has been used is the potential-field source 



surface (PFSS) model developed by Altschuler & Newkirk 
(1969), described below. 



Recently, ZDI maps for four classical T Tauri stars have 
shown that they possess complex fields that depart signifi- 
cantly from dipolar configurati ons (Donati et al. [12007 Do- 



nati et aL]|2008 Hussain et al.|[2009 ). This is significant as 
the field is thought to interact with circumstellar discs in 
ways that can depend on what kind of field the stars have. 
This interaction can have a large effect on the accretion of 
material in circumstellar discs onto stellar surfaces. Magne- 
tospheric accretion models for classical T Tauri stars usually 
assume that the star possesses a dipolar magnetic field ( |Ca-| 
menzind|1990[ [Koenigl 1991). More complex fields, however, 
can affect several aspects of accretion models. These effects 
include the locations of hotspots where accretion columns, 
funneled along field lines, impact onto the stellar surface 
( Gregory et al.|[2006 ) and angular momentum transfer be- 



tween discs and their stars (Long et al. 
et al.||2009t 



2008, Romanova 



Significant complexity is also typically found in the 
magnetic fields of main sequence stars that have been ob- 
served with this technique. One of the stars that has been 
most extensively studied in this way is the rapidly rotating. 



active star, AB Dor (e.g. Donati Collier Cameron||1997 



Donati et al.||1999[ [Jardine et al.||1999[ [Donati et al.||2003f . 



Jardine et al. ( 1999 ) and Hussain et al.| (|2000|) found that 
this technique used on observed radial field maps can accu- 
rately reproduce the observed azimuthal field maps up to a 



latitude of 60°. Jardine et al. (1999) suggested that the dis- 



agreement above 60° could be due to the presence of dark 
polar spots suppressing the Zeeman signature at these lat- 
itudes. They also noted that due to a tilt of about 60° to 
our line of sight in the star's rotation axis, a large part of 
the surface is out of view, and thus all studies on this star 
must be carried out without any information over about an 
entire third of the surface. 

While the ZDI technique has been applied with great 
success, allowing for the global structure of magnetic fields 
to be explored in stars other than the Sun, several funda- 
mental difficulties arise from the finite ability of the ZDI 
maps in picking up field structures. 

As stars are unresolved, differentiating between mag- 
netic features above and below the equator, when using the 
Doppler imaging and ZDI techniques, requires that the ro- 
tation axis of the stars be tilted to a certain extent. Unfor- 
tunately, this means that ZDI maps are unable to reproduce 
the field over the entire surface of any star. 

Another problem comes from the presence of dark re- 
gions on the surface of stars that are generally assumed to 
be the stellar analogue of sunspots. If this assumption is 
valid, then they would presumably contain strong magnetic 
fields with strengths that exceed the rest of the stellar sur- 
face. fSoIankT (2002) considers the validity of this assump- 
tion, perdyugina (2005) argues that the differences between 
observed starspot filling factors and magnetic filling factors, 
along with the fact that regions of intermediate brightness 
have been observed to contain the strongest magnetic fields 
on ZDI maps, indicates that the fields from starspot um- 



brae have yet to be detected. If this is the case, it may be 
significant limitation to ZDI maps and subsequently to the 
reconstruction of the 3-dimensional fields. 

The surface field maps also suffer from having limited 
spatial resolution. Since ZDI maps actually measure mag- 
netic flux averaged over each surface element, regions of 
opposite polarity will cancel each other out and thus not 
be picked up on the map if they are not resolved. The 
field strengths derived from Zeeman-Doppler Imaging are 
always lower than those derived from Zeeman broadening. 
Any structures not present on ZDI maps will obviously not 
be accounted for when the field is extrapolated to higher 
radii. 

The aim of this paper is to explore some of the lim- 
itations that these effects impose on extrapolating global 
coronal fields. Starting from the surface fields of two simu- 
lated stars, where these problems are not present, the global 
coronal fields and the emission measures are calculated. The 
effects discussed above are simulated so that the calculations 
can be repeated and compared. 



2 MODEL 

2.1 Simulated Stellar Magnetograms 

The radial component of the magnetic field at the base of 
the stellar corona is simulated using the magnetic flux trans- 
port model of Mackay et aL] (|2004|) and van Ballegooijen 



et al. ( 1998| ). The model describes the time evolution of the 
radial component of the magnetic field at the stellar surface 
through the magnetic induction equation under the effects 
of differential rotation, meridional flow, supergranular diffu- 
sion and the emergence of new magnetic fields in the form 
of magnetic bipoles. 

This paper will use two simulations of stars from 



Mackay et al. (2004). These stars differ by the range of lat- 



itudes at which flux is emerged and the meridional flow. In 
both cases, a comparable level of magnetic flux is emerged. 
The star given the greater meridional flow and range of lat- 
itudes for flux emergence (hereinafter the " active star" ) has 
flux emerging between latitudes of 10° and 70° in both hemi- 
spheres and a meridional flow of lOOms"^. The other star 
(hereinafter the "inactive star") has flux emerging between 
latitudes of 10° and 40° in both hemispheres and a merid- 
ional flow of llms~^. For both stars, the diffusion of mag- 
netic flux is characterized by the diffusion constant, D (see 
Equation 1 of [Mackay et al.^ 2004), which is given the value 
of 450/cm^s~^. The activity of the active star is analogous 
to the young active star AB Dor while that of the inactive 
star is analogous to the Sun. 

The magnetic activity of both stars is simulated to give 
periodic variations over 11-year cycles. Butterfly diagrams 
showing the latitudinal concentration of magnetic activity 
over these cycles for both stars are given in Fig. [l] As can 
be seen, the concentrations of polar flux visible in Fig. [2] are 
present throughout the cycles of both stars. Both stars show 
reversals in polar flux between every cycle. Sample radial 
field surface maps are given in Fig. |2] These correspond to 
times of maximum magnetic activity for both stars. Black 
contour lines are drawn at field strengths, \Br\^ of 500G. 
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Figure 1. Simulated butterfly diagrams for an active star (top) 
and an inactive star (bottom) showing the signed magnetic flux 
integrated over all longitudes as a function of latitude and time. 
Both diagrams show almost two and a half 11-year activity cycles. 
Red and blue regions represent net fluxes of opposite polarity. 



2.2 Potential-Field Extrapolation 

The simulated surface fields are extrapolated to higher radii 
using the potential-field source surface model described by 



Altschuler & Newkirk (1969). The field is assumed to be 



current-free (V x B = 0) and thus satisfies Laplace's equa- 
tion (V^$ = 0). The 3D coronal field can then be described 
as the sum over multiple spherical harmonic components 
and is given by the following equations where Pim are the 
associated Legendre polynomials. 
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Figure 2. Surface radial fleld maps for an active star (top) and 
an inactive star (bottom) at a time of a relative maximum in 
magnetic activity. Black lines give contours at 500G. Red and 
blue regions represent magnetic flelds of opposite polarity. 



at a radius of 3.4i?*. With this, the radial component of the 
field becomes 



mcf) 



(1) 



where 



fir J) 



-(/ + l)(^)-a + 2)^^(^)-(2. + l)(^y-r 
/ + l + /(^)-(2^ + l) 



The surface radial field maps are used as a lower bound- 
ary condition to derive the values of cim for all values of / 
and m under consideration. For this study, values of / up 
to 31 will be considered. The field is extrapolated using a 
modified version of the code developed by |van Ballegooijen| 
eTaTl ([19981. 



Once the global field is determined, the coronal gas pres- 
sures can be estimated using the model described in |Jardine| 
et al. ( ^002 ) and Gregory et al. ( 2006 ). The gas pressure in 
regions corresponding to open field lines is taken to be zero 
and thus only regions of the stellar corona confined by closed 
magnetic loops contribute to the emission measure. The gas 
pressure at the stellar surface is assumed to be proportional 
to the magnetic pressure (po oc Bq) and to vary along field 
lines assuming the gas is in hydrostatic equilibrium, given 
by 



p — poexp 



m 
kf 



Qsds 



As an upper boundary condition, we assume that the 
field becomes radial at a source surface which is taken to be 



where Qs is the component of effective gravitational accelera- 
tion along the field line. Once a temperature is assumed, the 
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density rie can then be determined. We assume an emissiv- 
ity proportional to rie to obtain the X-ray emission measure. 
We consider both its magnitude and rotational modulation 
throughout the stellar magnetic cycle. 

2.3 The effects of star spots, the inclination of the 
stellar rotation axis and finite surface 
resolution 

Once the coronal magnetic field structure, gas pressures and 
emission measures have been calculated using the full surface 
magnetograms, the process is repeated using magnetograms 
that have been degraded to mimic the effects of a) the mask- 
ing of the Zeeman signature in dark starspots that leads to 
a reduced magnetic flux measurement b) the lack of infor- 
mation from a hemisphere that is out of view because of the 
inclination of the stellar rotation axis and c) a finite surface 
resolution. 

The presence of star spots is simulated by cutting from 
the maps all regions with fields, \Br\^ stronger than a thresh- 
old of 500G. Similarly, the obscuration of a portion of the 
star's surface due to a tilt in its rotation axis is simulated 
by cutting out an entire hemisphere. This corresponds to a 
limiting case where the star's rotation axis is parallel with 
our line of sight and we can only see one hemisphere. If 
a star actually had this tilt angle, the features on its sur- 
face would not be Doppler shifted as the star rotates, and 
thus stellar magnetograms would be impossible to produce. 
Even though this does not correspond to a physically real- 
istic situation, we nonetheless consider it as it provides the 
most severe limitation on the fraction of the star that is vis- 
ible. We note however, that while the very early ZDI maps 
were produced by reconstructing the radial component, or 
all three vector components, of the field at each point on 



Star 


A^spots (%) 


A*^hemisphere 


A^smeared (%) 


Active 


19 


52 


17 


Inactive 


11 


51 


10 



the surface independently ( Brown et al.||1991 Donati et al. 
1999), the more recent reconstructions ( [Donati et al.||2001 



Donati et aT1|2QQ6[ [Morin et ar]|2008t describe the field in 
terms of spherical harmonics whose amplitudes are fitted to 
the line profiles. This allows the reconstruction of magnetic 
modes that are symmetric or antisymmetric about the equa- 
tor and therefore some indication of the global symmetry of 
the field can be obtained by comparing the fit to the data 
of the modes, even if some portion of one hemisphere is out 
of view. 

The opposite limiting case, where the star's axis of ro- 
tation is perpendicular to our line of sight, also makes it im- 
possible for ZDI maps to be produced. In this case, features 
at the same latitude in opposite hemispheres are indistin- 
guishable in the broadened spectra and so it is impossible 
to assign magnetic elements to a single hemisphere. 

The effect of having low-resolution ZDI maps is sim- 
ulated by smearing the surface field such that the original 
maps, with a latitudinal resolution of just over 1° and an 
equal longitudinal resolution at the equator, are replaced by 
maps with latitudinal resolution of around 11° degrees and 
a longitudinal resolution of around 8° at the equator. As 
ZDI maps measure magnetic flux, each new surface element 
contains a field which is simply the area weighted mean of 
the fields in the old elements contained within it. 

In all three of these cases, the magnetic field is ex- 
trapolated and the results compared to the original case. 
In analysing the difference between the results, two aspects 



Table 1. Percentage drop, averaged over one activity cycle, in 
total magnetic flux, <I>, over the stellar surface for both stars. This 
drop is due to the suppression of magnetic field in dark starspots 
and an entire hemisphere out of view and the cancellation of 
opposite polarity magnetic fields due to the smearing of the stellar 
surface magnetogram. 



will be considered. These are the effects on the structure of 
the magnetic field and the effects on the coronal gas pres- 
sures, densities and emission measures. Also, the effects on 
the open flux and the field lines that can potentially support 
accretion from a circumstellar disc are explored. 



3 RESULTS 

3.1 The surface maps 

The effect of degrading the surface magnetograms can be 
seen very clearly by plotting the amplitudes cim of the spher- 
ical harmonic components from equation ([T]) as a function 
of / and m. This is shown in Fig. [3] for a period of maxi- 
mum activity during the magnetic cycle. This corresponds 
to the same time as that shown in the magnetograms in Fig. 
|2] The first column shows the distribution of power across 
the modes for both the active star (top) and the inactive 
star (bottom). The magnetic field of the active star shows 
power across a much greater range of modes, corresponding 
to its more complex and widely-distributed field. The sec- 
ond column, illustrating the effect of dark star spots, shows 
the influence on the field structure of cutting out all regions 
of field above 500G. This clearly has a greater effect on the 
field of the more active star since it has a greater coverage 
of these strong- field regions. While some power is removed 
from all the modes, it is the high-order modes that are most 
strongly suppressed, while the dominant low-order modes 
are largely unaffected. The third column illustrates the ef- 
fect of removing field from the unobserved hemisphere. This 
has a significant effect on the field structures of both stars. 
It allows power to be distributed more evenly over the differ- 
ent modes, giving the effect of blurring the Qm plot given in 
the second column of Fig. |3] The fourth column illustrates 
the effect of smearing the original magnetogram. This re- 
moves some of the higher order components while leaving 
the lower order components unchanged. As the strength of 
higher order components decreases with increasing radius 
faster than lower order components, this affects the field on 
smaller scales more than on larger scales. 

Table [l] shows the percentage drop (over the stellar ac- 
tivity cycle) in the total unsigned surface magnetic flux for 
both the active and inactive star. The case where an entire 
hemisphere is cut from the surface magnetograms presents 
the largest drop in magnetic flux of about 50% on average. 
For the active star, the case where the surface magnetogram 
is smeared gives a large drop in magnetic flux due to can- 
cellation of opposite polarity flux. 

Although these stars show large unsigned magnetic flux 
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Figure 3. Plots showing the strength of each cim coefficients given in Equation^ Results are shown for the full data case (first column), 
the field cutoff above 500G (second column), the field cutoff in an entire hemisphere (third column) and the field smeared (forth column) 
for both the active star (top) and inactive star (bottom). 



over the entire magnet ogr am, according to Gauss' law for 
the magnetic field (^^ B • = 0) the total signed flux 
must always be zero (in the absence of magnetic monopoles). 
However, when regions of the magnetograms are suppressed, 
if more flux of one polarity is lost, this condition will not 
be satisfied. This is the case when the fields stronger than 
500G and when the field in an entire hemisphere are sup- 
pressed, although in both cases, the total signed flux is usu- 
ally smaller than 2% of the total unsigned flux. In these cal- 
culations, the surface magnetograms are represented as the 
sum over many spherical-harmonic components. The total 
signed flux for each individual component is zero and thus 
for the spherical-harmonic expansion of the surface magne- 
tograms is also zero. Thus, the extrapolated magnetic fields, 
in most cases, cannot perfectly represent the original surface 
magnetograms when regions of the magnetograms are sup- 
pressed. Instead of being zero, the total unsigned flux in 
the hemisphere that is suppressed after the extrapolation, is 
usually around 5 - 10 % of its original (unsuppressed) value. 

3.2 The magnetic field structure, accretion and 
open fiux 

Extrapolated magnetic fields for both stars are shown in 
Fig. [4] These dia grams give the field structure at a time of 
maximum activity corresponding to the surface maps given 
in Fig. [2] 

One of the reasons why understanding the structure 
of stellar magnetic fields is important is that it will help 
to understand magnetospheric accretion to young pre-main- 
sequence stars from their circumstellar discs. In this section, 
we will discuss the consequences of missing information in 
this paper has on calculated accretion filling factors and the 
latitudes at which accreting field lines attach to the stellar 
surface. For both stars and at all times, the effect of low 
resolution surface magnetograms is insignificant and will not 
be discussed here. 

The locations on the stellar surface where accretion can 
take place are calculated by tracing field lines from every 
grid point on the stellar surface. For these calculations, the 
global magnetic fields are extrapolated assuming a source 
radius of lOi?*. Only the field lines that cross the equator 
and extend to radii between 0.5i?co and l.Oi^co, where Rco 



is the equatorial corotation radius, are assumed to be able 
to support accretion. All of these field lines are assumed to 
be supporting accretion. In this case, Rco (given by Rco = 
{GM^/u'^Y^^) is calculated assuming the star has solar mass 
and radius and a rotation period of 2.5 days and has a value 
of approximately 8i?* . 

Fig. |5] gives a butterfly diagram showing the latitudes 
at which field lines that can support accretion connect with 
the stellar surface over most of an activity cycle. When the 
field with strengths above 500G are suppressed, it can be 
seen that there is only a small effect on the accretion filling 
factors for the inactive star. The active star shows small in- 
creases at some times during the cycle with an increase in 
the amount of low latitude regions contributing to accretion 
at most times in the cycle. The most significant effect on the 
accretion filling factors and latitudes comes in the case where 
the field is suppressed in an entire hemisphere. This is clear 
for both stars and is due to the fact that without regions of 
stronger and weaker field in the invisible hemisphere, field 
lines that are being traced from the visible hemisphere can 
connect to anywhere with no preferred location. However, 
when an entire hemisphere is out of view, the accretion pat- 
tern in the visible hemisphere tends to remain the same. 

Another important aspect that the magnetic structure 
can have on accretion is its effect on the magnetic pressure 
in the equatorial plane, which determines the fields ability 
to disrupt the circumstellar disc. Fortunately, it has been 
found that for both stars, the missing information causes no 
significant change in the strength of the magnetic field in the 
equatorial plane, and its variation as a function of distance 
from the surface. 

Finally, the effect that the missing information has on 
calculated open magnetic flux will be explored briefly. The 
two important aspects that will be discussed are the fraction 
of open flux to total flux, and the latitudinal distribution. 
These two aspects can affect the rate at which stars lose 
angular momentum through stellar winds. 

In every case, the missing information causes a large 
drop in the total open flux due to the drop in total unsigned 
flux. Table [2] gives the average percentage drop in the frac- 
tion of open flux to total flux for each case. It can be seen 
that in almost every case the missing information causes an 
increase in the fraction of open flux to total flux. The effect 
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Figure 4. Extrapolated coronal fields for the active star (right) and the inactive star (left). Presented for the full data (first row), the 
field cut off above 500G (second row), the field cutoff in an entire hemisphere (third row) and the field smeared (fourth row). Open field 
lines are given in blue and closed field lines in red. These images correspond to extrapolations from the surface field maps given in Fig. 
I2]and thus show the star at a period of maximum magnetic activity. 
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Figure 5. Butterfly diagrams for the active star (top row) and inactive star (bottom row) for the fuh case (left column), suppression 
of the field stronger than 500G (middle column) and suppression of the field in an entire hemisphere (right column). The black regions 
show the latitudinal locations where field lines able to support accretion from a circumstellar disc attach to the stellar surface. 
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Active 


15 
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15 
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Inactive 


9 


-2 


20 
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-68 


-1 
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Table 2. Percentage increase, averaged over one activity cycle, 
in the fraction of open fiux to total fiux, / = /\{(^ open total) ^ 
for fields with missing information due to starspots, an entire 
hemisphere out of view and a limited resolution for both active 
and inactive stars. 



Table 3. Percentage increase, averaged over one activity cycle, in 
coronal electron density, rie, averaged over the volume of plasma 
in closed magnetic loops. Results given for both stars due to the 
absence of magnetic fields in starspots, a hemisphere out of view 
and the effects of having a surface magnetogram with a reduced 
resolution. 



of the missing information not only decreases the open flux 
and the total flux, but also changes the nature of that flux. 
Typically the flux that is lost is closed. The loss of informa- 
tion leads to a greater fraction of the flux through the stellar 
surface being open. 

For both the active and inactive stars, the open flux 
tends to be located at the stars' poles. The latitudinal dis- 
tribution of the open flux is not generally affected by the 
missing information for both stars. An exception to this is 
the case of the suppression of the field stronger than 500G 
on the active star. In all other cases, the open flux is gener- 
ally located around the poles, whereas in this case, the open 
flux is generally located at slightly lower latitudes. 



3.3 The coronal density and emission measure 



As described in section [2^2] the density at each point in the 
corona is determined assuming isothermal, hydrostatic bal- 
ance along the field line through that point. The one free 
parameter in this process is the pressure at the base of 
each field line. Since this is assumed to be proportional to 
the magnetic pressure there, any degrading of the surface 
magnetic map will affect the modelled coronal pressure and 
hence the modeled density. This can be seen clearly in Fig. 
|6] which shows histograms of the number of volume elements 
in the corona at specific electron densities for a time close 
to cycle maximum (corresponding to Fig. [5] and Fig.js]). The 



results are presented for both stars and all cases under con- 
sideration. The most significant impact on the distribution 
of pressures and hence densities in the corona is for the case 
where the spot fields are removed. Since this cuts out the 
regions of strongest field, it also cuts out the highest den- 
sity values, as can be seen by the loss of the high-density 
tail in the second column of Fig. [6] This would be detected 
observationally as an enhancement of the densities derived 
from X-ray lines (e.g. Testa et al. 112004 ) compared to the 
model values (for the same emission measures). The frac- 
tion by which the models would underestimate the densities 
as a result of the loss of information in the magnetograms 
is shown in Tabled It can be seen that in the cases of miss- 
ing information due to starspots and a limited resolution, 
large decreases in the average electron densities are present. 
In the case of missing information in an entire hemisphere, 
there is no large decrease in the electron densities with the 
inactive star remaining constant and the active star showing 
a significant increase. 

This alteration of the modeled coronal density is a re- 
sponse to changes in the surface magnetograms and is also 
apparent when the emission measures are calculated. Fig. 
IXl shows the variations of emission measures and the rota- 
tional modulation of the emission measures for just over one 
11-year activity cycle for both stars. In all cases, the vari- 
ation in emission measure closely follows the variation in 
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Figure 6. Histograms showing the number of regions in the corona with specific electron densities for the active star (top row) and the 
inactive star (bottom row). The cases of the full data case (first column), the field cutoff above 500G (second column), the field cutoff 
in an entire hemisphere (third column) and the field smeared (forth column) are given. 
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Figure 7. Emission measures over one 11-year activity cycle for the active star (top row) and inactive star (bottom row). Results are 
shown for the full data case (first column), the field cutoff above 500G (second column), the field cutoff in an entire hemisphere (third 
column) and the field smeared (forth column). The thickness of the lines show the rotational modulation with the tops of the lines 
corresponding to the maximum and the bottoms of the lines corresponding to the minimum in visible emission measures as the stars 
rotate. 



total magnetic flux over the activity cycle shown in Table 
[l] For the active star, the variation of the emission measure 
is fairly robust and is only significantly affected by the re- 
moval of the field in the spots shown in the second column. 
The very large drop in the emission measure is due to the 
large fraction of the surface covered in fields greater than 
500G. The inactive star has many fewer high-field regions 
and consequently a much larger rotational modulation. It is 
also much more severely affected by both the removal of the 
spot fields and also by the removal of the field in the lower 
hemisphere. In both these cases, the form of the variation is 
changed (indeed the cyclic nature of the variation is all but 
lost when the spot fields are removed). For both the active 
and inactive star, the smearing of the maps has least effect, 
while the removal of the spot fields has the greatest effect. It 
is only for the inactive star, however, the the form of the cy- 
cle is changed. This suggests that calculating stellar coronal 



emission measures from ZDI maps may fail to reproduce the 
cyclic variations in emission measures for some stars even if 
magnetic cycles are present. 

The large drops in the modeled emission measures due 
to the suppression of the strongest fields from the surface 
magnetograms clearly depends on the spot filling factors. 
This is shown in Fig. [8] for both stars. It can be seen that 
the dependance on spot filling factors is the same for both 
stars. 



4 CONCLUSION 

A potential field source surface model has been used to ex- 
trapolate the global coronal magnetic fields from surface ra- 
dial field maps of two simulated stars. Using these extrap- 
olated global fields, the coronal emission measures and ro- 
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in starspot umbrae (Berdyugina 2002), Berdyugina et al 



(2008) was able to detect Stokes V signals on at least four 
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Figure 8. Drop in calculated emission measure for both active 
and inactive stars against the spot filling factor (fraction of stellar 
surface area with field above 500G) for the active star (green) and 
inactive star (red). Separate data points correspond to different 
times over the stellar activity cycle. The thin green line is a best 
fit to the data for the active star and is given by log(A£^M) = 
53.95/g •^■'^ — 0.11. The thin red line is for the inactive star and is 
given by \og{/\EM) = 53.42/0 01 _ 0.53. 



tational modulations were calculated over ll-year activity 
cycles. The calculations were repeated for several cases to 
simulate the absence of information due to starspots, lim- 
ited spatial resolution and limited surface visibility due to a 
tilt in a star's rotation axis. The results were compared in 
order to understand the limitations of reconstructing coro- 
nal magnetic fields in the absence of this information. In 
all cases, the field structures and emission measures were 
aflFected to various degrees. 

In the presence of the simulated starspots, the global 
flux levels calculated from the surface magnetograms are 
reduced, but the fleld structure is not drastically aflFected 
on large scales. The largest change in the fleld comes on a 
smaller scale from the loss of the most active regions. These 
results suggest that the presence of starspots does not signifl- 
cantly aflFect the extrapolation of large scale coronal fleld but 
is relevant for studying smaller scale fleld structures. This 
does however lead to the largest decrease in the emission 
measure from the original case of the three effects studied, 
hinting that starspots may be one of the largest problems 
when reproducing stellar emission measures. The rotational 
modulation of emission measures are also signiflcantly af- 
fected when the strongest fleld regions are lost. 

The loss of fleld measurements in dark spotted regions 
has the greatest impact on the gas pressures inferred from 
the magnetic fleld extrapolation. In particular, it reduces 
the maximum gas pressures (and hence densities) that can 
be predicted. 



Using the doppler imaging technique (Vogt & Penrod 



1983), it is possible to determine the positions and sizes of 



starspots of suflficient sizes on the surfaces of stars. One pos- 
sible solution to the problems that starspots introduce may 
be to artiflcially impose fleld onto the surface magnetograms. 
Another possible solution may be to detect the flelds within 
starspots directly. Using molecular lines that are only formed 



M dwarfs. 

The missing information due to a lack of visibility of 
an entire hemisphere has a greater effect on the flnal fleld 
structure than that of the loss of high fleld strength regions. 
In this case, it appeares that the fleld is weaker and more 
spread out over different spherical harmonic components. 
The estimated coronal gas pressures and emission measures 
are not signiflcantly affected thought in the stellar activity 
cycle and on average half of the high pressure regions of the 
stellar corona are preserved. However, the emission measures 
drop signiflcantly on the rare occasions during the stellar 
cycle when most of the magnetic flux happens to be in the 
hemisphere out of view. This is a problem for estimating the 
emission measure from real stars as one will not be able to 
know how much magnetic flux covers the missing regions of 
the observed star. 

The observed drop in emission measure for both these 
cases is due to a loss of the highest coronal gas pressure 
regions which flU the volume above the high fleld strength 
regions. The modeled gas pressures depend on the assump- 
tion that the pressure at the base of closed fleld lines is 
proportional to the magnetic pressure {po oc Bq). Due to 
this assumption the variations in gas pressures and emission 
measures for both stars follow closely the total magnetic 
flux. If a weaker dependance on the magnetic pressure is 
taken {po oc Bq where a < 2), then the dependance of the 
emission measure on the regions of high fleld strength will 
be reduced. Of course a stronger dependance {a > 2) would 
have the opposite effect. 

The flnal case of taking into account the flnite resolu- 
tion limitation on observational reconstructions of surface 
magnetic flelds was shown to have a simplifying effect on 
the extrapolated fleld. This suggests that the resolution is 
only an issue when studying highly complex flelds. 

It can be seen in section [3^ that calculations involving 
fleld lines that cross the equator and extend to a radius that 
could coincide with the inner region of a circumstellar disc 
can be affected in various ways. When the fleld in an entire 
hemisphere is suppressed, the modeled locations of hotspots 
are largely unaffected in the opposite hemisphere. We are 
unable to reproduce the location of hotspots in the hemi- 
sphere that is not visible leading to an increase in the calcu- 
lated surface flUing factors of hotspots in that hemisphere. 
In order to study accretion onto stars using fleld structures 
extrapolated from ZDI maps, these results suggest that it 
is better to reproduce plausible fleld conflgurations in the 
hemisphere out of view (as in Gregor y et al.|20 06 ) although 
it is possible to reproduce the location and distribution of 
accretion hotspots in the visible hemisphere, even with in- 
complete surface magnetograms from ZDI maps. 
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